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ABSTRACT 

We present a new method that deals with the uncertainty in matter-clustering in 
cosmic shear power spectrum analysis that arises mainly due to poorly understood 
nonlinear baryonic processes on small-scales. We show that the majority of information 
about dark energy physics contained in the shear power comes from these small- 
scales; removing these nonlinear scales from a cosmic shear analysis results in a 50% 
cut in the accuracy of measurements of dark energy parameters, marginalizing over 
all other parameters. In this paper we propose a method to recover the information 
on small-scales by allowing cosmic shear surveys to measure the nonlinear matter 
power spectrum themselves and marginalize over all possible power spectra using 
path integrals. Information is still recoverable in these nonlinear regimes from the 
geometric part of weak lensing. In this self-calibration regime we find we recover 90% 
of the information on dark energy. Including an informative prior, we find the nonlinear 
matter power spectrum needs to be accurately known to 1% down to fc = 50/iMpc~^ 
to recover 99% of the dark energy information. This presents a significant theoretical 
challenge to understand baryonic effects on the scale of galaxy haloes. However self- 
calibration from weak lensing may also provide observational input to help constrain 
baryon physics. 
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1 INTRODUCTION 



Cosmic shear has been identified as being a particularly sen- 
sitive tool in understanding dark energy (Albrecht et al., 
2006; Peacock et al., 2006), dark matter (Massey, Kitching, 
, Richards, 2010), neutrino physics (Hannestaad, 2010) and 
■ potential depatures from general relativity. There are a num- 
ber of on-going (CFHTLenS, Pan-STARRS) and planned 
experiments (KIDS, DES, LSST, Euclid) whose primary sci- 
entific goals are to use cosmic shear to constrain cosmologi- 
cal parameters. 

However, as we show in this article, the majority of 
dark energy information comes from small scales, that are 
expected to be influenced by poorly understood non-linear 
effects such as baryonic processes (see White, 2004; Zhan & 
Knox, 2004, Mead et al., 2010; Rudd et al., 2008; Guillet et 
al., 2010 for example); which may be very difficult to model 
to sufficient accuracy. Indeed we may not even know the 
sign of contribution from baryonic effects. If dissipation is 
the main effect, we can imagine baryonic collapse will lead 
to an enhancement of matter clustering, or Baryonic Com- 
pression. But nonlinear feedback from star or AGN produc- 
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tion could also blow out baryons, reducing the overall mass 
and allowing the dark matter to disperse. The formation of 
stars or AGN in simulations is usually governed by sub-grid 
semi-analytic prescriptions, which makes the understanding 
of such processes uncertain. As a step forward it would be 
very useful to even quantify our uncertainty on the physics 
involved on different scales. 

In addition to the uncertain physical effects there is the 
practical consideration of the accuracy of fitting formulae. 
The matter power spectrum as a function of redshift can be 
computed using linear perturbation theory of the underlying 
initial density field, but by comparison with N-body simula- 
tions it is apparent that at scales where the matter overden- 
sity 5 becomes greater than 0.1 the linear theory predictions 
cannot be used as non-linear effects in the growth of struc- 
ture become dominant. The most widely used corrections are 
Smith et al. (2003) halofit. Peacock & Dodds (1996) and 
the Coyote formula Heitmann, et al, (2010); halofit and 
Coyote are accurate to approximately 5 — 10%. In addition 
to this uncertainty the formula are proposed by only sam- 
pling a small discrete number of points in parameter space, 
and currently do not include baryons. 

Finally, to extract cosmological information from the 
statistical properties of the density field means we require a 
covariance matrix. Nonlinear growth of dark matter struc- 



© 2010 RAS 



2 T. D. Kitching & A. N. Taylor 



ture will already correlate estimates of the matter power 
spectrum on different scales (e.g. Kiessling, et al. 2010), and 
including baryonic physics we can expect even stronger co- 
variances between scales due to feedback processes. 

The cosmic shear power spectrum depends on the mat- 
ter power spectrum through an integral over the line of 
sight distance, with a geometric lensing kernel. The lensing 
power spectrum contains cosmological information, through 
the lensing kernel, even on small scales. Rather than throw 
this information away, wc propose that one can include oven 
very uncertain non-linear scales in a lensing analsysis, if one 
correctly marginalizes over the uncertainty. Using the path- 
integral marginalization techniques, presented in Taylor & 
Kitching (2010) and Kitching & Taylor (2010), we derive 
an expression, in the self-calibration regime, for the cosmic 
shear covariance that includes the geometric lensing kernel 
information from all scales. 

If an informative upper bound on the functional be- 
haviour of the non-linear power spectrum can be determined 
(from simulations for example) the residual uncertainty in 
the functional form of the power spectrum can be marginal- 
ized over simultaneously with the cosmological parameters 
of interest. We place a requirement on the external fuctional 
prior such that the cosmological constraints from future all- 
sky cosmic shear experiments are not degraded below a level 
needed to determine dark energy physics. 

This article is organised as follows, in Section 2 we 
present the marginalized likelihoods and Fisher matrices for 
the non-linear tomographic cosmic shear power spectrum. 
In Section 3 we present results, and in Section 4 we present 
our conclusions. 



2 METHOD 

2.1 Cosmic Shear Tomography 

In this article wc will focus on tomographic cosmic shear, 
where the shape information of galaxies is used, and objects 
are binned by their estimated redshift, to create a series 
of 2D cosmic shear maps. The auto- and cross-correlation of 
these maps can be used to generate a series of power spectra 
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We follow the notation of Joachimi & Bridle (2009). The 
comoving distance is r, rn is the horizon distance, while 
'S'fe(r) = sin(r), r, sinh(r) for curvatures k — — 1,0,-|-1, a is 
the scale factor and Pss{£/Sk{r);r) is the 3D density matter 
power spectrum. The comoving galaxy probability distribu- 
tion is given by Piir). The ij subscripts refer to redshift 
bins, where the shear field is approximated as a series of 
correlated 2D planes. This can be generalized to a full 3D 
shear estimator (Kitching, Heavens, Miller, 2010). We will 



neglect intrinsic alignment terms so that the covariance of 
the cosmic shear can be written as 

c^.(£) = 2Cim{e)Cii{e) (4) 

where fi = ij and v = ml refer to redshift-bin pairs. 

2.2 Information in the Non-Linear Regime 

One subtlety, a result of the Limber approximation used in 
equation (1), and explicitly highlighted by Kitching, Heav- 
ens & Miller (2010), is that the radial fc- modes in the mat- 
ter power spectrum are associated with the azimuthal i- 
modes through I — kr{z) where r{z) is the comoving dis- 
tance for redshift z. Commonly angular wavenumbers of 
up to ^max = 20,000 are used in theoretical predictions of 
weak lensing power. If no cut in Fourier 3D wavenumber 
k is applied for large I modes the resulting scales probed 
in the matter power spectrum can extend into the very 
highly non-linear regime. For instance for ^max = 20,000 
and r{z = 0.1) = 400Mpc//i we find k = 70/iMpc~\ or a 
physical scale of 126Kpc//i, within the dark matter halo of 
gala:xies. We demonstrate this in Figure 1 where we show 
the (^, fc) plane and the accessible modes, that lie on the 
lines for a 15-bin tomography experiment from redshift 0.1 
to 1.5; it is clear that all modes with I > 2000 probe the 
non-linear regime exclusively. 

In this article we will make defined cut in radial fe-mode, 
where k < fcmax for all tomographic bins. This allows scales 
to be probed in a controlled manner so that the non-linear or 
very non-linear regimes can be removed from the analyses. 
Equivalently, this can be intepreted as a redshift-dependent 

I cut where £max(2) = fcmaxr(«). 

In Figure 2 we show the predicted dark energy Figure of 

Merit ^ , calculated using a Fisher matrix (Tegmark, Taylor & 
Heavens, 1997; see Hu 1999 for the cosmic shear tomography 
Fisher matrix), as a function of the maximum fc-mode used 
in the matter power spectrum for a Euclid-like^ tomographic 
cosmic shear survey, and a set of cosmological parameters 
that includes curvature"'. It is clear that the majority of the 
dark energy information comes from modes in the non-linear 
regime (although this statement contains some uncertainty 
because we have used the Smith et al, 2003 halof it correc- 
tion). 

We find convergence at fcmax ~ 50ftMpc~^. This is in 
agreement with previous studies (e.g. Dore, Tingting, Ue-Li, 
2009). We can understand the behavior of this plot where 
the increase in FoM around k = 0.1/iMpc~^ corresponds to 



^ FoM; Albrecht et al., (2006), defined as the area con- 
strained in the uncorrelated (wQ,Wa) plane with FoM = 

-^/['\/(^w0u)0^'"O""» ~ (-'^■uiOwa) )]• 

^ Refregier et al., (2010), 20,000 square degrees, with median 
redshift of Zm = 1.0, 35 galaxies per square axcminute, a photo- 
metric redshift uncertainty of 0.03(1 + z). We assume the galaxy 
number density is given by n(z) oc z'^cxp{—lAz/zmY"^, and use 
10 redshift bins. Throughout we include no priors. 
^ Throughout wc will use a cosmological parameter set that al- 
lows for curved cosmologies with parameters f2m, Vl^^, f!;,, h, 
o"8, and ris given by (0.25,0.7,0.75,0.8,0.95), and parameterize 
the dark energy equation-of-state using a first-order Taylor ex- 
pansion, w{z) = two + (1 — a)wa (Linder, 2003; Chevallier and 
Polarski, 2001) with {wo,Wa) = (-0.95,0.0). 
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Figure 1. The k modes that are in principle accessible to a 
tomographic experiment. As an example we use a 15 bin tomog- 
raphy experiment with bins at O.ln from z = 0.1 to 1.5; solid lines 
drawn are £ = kr{z). Tomography uses modes that lie on the lines 
shown. 3D cosmic shear (Kitching, Heavens, Miller, 2010) would 
use all modes below the green line, which marks the maximum 
redshift of the example experiment. The dashed line shows the 
approximate non-linear cut off of fc 0.5/i.Mpc~^. 
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Figure 2. The dark energy Figure of Merit (FoM) as a func- 
tion of the maximum fc-mode used in the tomographic cosmic 
shear Fisher matrix analysis for a Euclid-like survey, with no 
external priors included. The maximum £ range in the C{£) 
power spectrum is calculated as a function of redshift using 

^max{2) = min[fcmaxr'(2), 2 X lO^]. 



the nonlinear regime for the matter power-spectrum. The 
flattening between k = 0.5 and k — 3/iMpc~^ occurs near 
the linear peak of the matter power spectrum which is less 
sensitive to dark energy. A second increase in information 
appears around the peak of the nonlinear regime, and then 
again saturates at very high k. 



2.3 Path Integral Marginalization 

The main issue we are concerned with in this paper is the 
uncertainty in the shape of the nonlinear matter power spec- 
trum due to baryon effects. In the absence of an accurate 
estimate of the nonlinear matter power spectrum one can ei- 



ther attempt to self-calibrate it from the cosmic shear data 
itself, or apply an informative external prior on its shape 
from hydrodynamical simulations. The self-calibration ap- 
proach finds the maximum likelihood solution to the free 
functional nonlinear matter power spectrum using the shear 
data alone. For a Gaussian likelihood for the matter power 
spectrum this can be done in a single-step using Newtons 
Method or a quadratic estimator (Taylor & Kitching 2010), 
or iteratively. We then analytically marginalize over all func- 
tional forms for the power spectrum permitted by the data, 
again assuming Gaussian-distributed nonlinear power. The 
resulting marginalized likelihood function for cosmological 
parameters is now independent of any initial fiducial choice 
of the nonlinear matter power, since we have found its max- 
imum likelihood value. This method can be applied to an 
arbitrary likelihood function for the data, but if we have 
Gaussian-distributed data, the resulting marginalized likeli- 
hood function for the data is again a Gaussian with a new 
covariance matrix. 



where C is the orginal covariance and 
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for functions ^pa{x), and we assume the functional Fisher 
matrix, Fap{x,x') for the function ipa{x), is invertible. fj, is 
the mean signal, which in our case is the dj {(.) tomographic 
power spectrum. 

If we assume an informative external prior, for exam- 
ple, if we know something about how baryons will affect the 
matter power spectrum from hydrodynamical simulations, 
we can include this information by constraining its shape to 
some accuracy which in general will depend on scale. This 
additional information can be included in our formalism by 
multiplying the likelihood with a Gaussian prior with some 
covariance Capix' ,x"). After marginalization, we again find 
an analytic expression for a general likelihood. In the case of 
Gaussian-distributed data this also modifies the marginal- 
ized covariance, C^^. This in turn can be simplified using 
the Woodbury relation to find that the marginalized data 
covariance is the original covariance C with a new term 
added. 



Cm = C + 



Slljaix') 5tpfl{x") 



(7) 



In the following subsections we derive an expression for 
the covariance of the tomographic cosmic shear power spec- 
trum in each case. Given the modified covariance C*^ a like- 
lihood function can then be constructed that accounts for 
unknown functional behaviour in the matter power spec- 
trum 



£o = AD^[C 
where AD„ = C. 
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data analysis. We assume a Gaussian likelihood, but in fact 
any likelihood estimator can be used (see Taylor & Kitching, 
2010). For predictive forecasts an associated Fisher matrix 
can be constructed assuming that the cosmological infor- 
mation is the mean tomographic power spectrum at each 
redshift 
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for a set of cosmological parameters, 0. 
2.4 Self-Calibration 

In the self-calibration regime the cosmic shear data itself is 
used to measure the non-linear power spectrum simultane- 
ously with cosmological parameters. To calculate the impact 
of this self-calibration we first find the functional derivative 
of the lensing tomographic power spectrum with respect to 
the matter power sepctrum, given by 

■n 



— ^*ij,j>iV) f 



(10) 



SPs6{i'/r,r) 

For a flat-prior in function-space over the non-linear mat- 
ter power spectrum, where we use the data itself to fit and 
marginalize over uncertainty in the matter power spectrum, 
the marginalized functional covariance is given by 

Cf^Ai) = [C^X^) - C-l, {£)P,,., WCMl)]-' (11) 
where 



Jo 



drdr' W^°{r)[F{£; r, r'T^W^'' {r'), (12) 



and the functional Fisher matrix for the matter power spec- 
trum is 



F(£;r,r') = ^ C;J(£)W^^°'^(r)W^°°(r'). 



(13) 



In practice this functional Fisher matrix is binned in redshift 
and scale using 200 bins in each dimension in this paper. 

After marginalizing over the nonlinear matter power 
spectrum, there is still cosmological information in the lens- 
ing power spectrum through the lensing kernel. The func- 
tional Fisher matrix in equation (13) picks out the geomet- 
ric dependency of the tomographic cosmic shear power spec- 
trum. Marginalizing over the functional matter power spec- 
trum has some similarity to the geometric shear-ratio test 
(Jain & Taylor 2003, Taylor et al 2007, Kitching et al 2007), 
where the shear signal behind galaxy clusters and groups is 
ratioed, and the mass drops out leaving just the geometric 
part of the lensing signal. This method can be used in the 
fully nonlinear clustering regime as there is no sensitivity to 
clustering. 

2.5 External Prior 

For an external prior we take the functional derivative of the 
tomographic cosmic shear power spectrum with respect to 
the matter power spectrum (equation 10) and, assumming 
a Gaussian distribution in function-space can now write a 
covariance for the tomographic cosmic shear power spectrum 
that marginalizes over all unknown functional behaviour of 
the matter power spectrum 

C^M) = C^^i) + / dralWr, r)W,'^°{r)Wit''{r), (14) 
Jo 

where a%{£/r,r) = {\SPss{£/r,r)\'^) is the functional scatter 
in Pss and we have assumed the functional covariance is 
diagonal in £ and r. 



3 RESULTS 

Here we present cosmic shear Fisher matrix forecasts, 
marginalizing over uncertainty in the matter power spec- 
trum, for the self-calibration and external prior cases de- 
scribed in Section 2. 



3.1 Self-Calibration 

In Figure 2 we show the FoM for the case that the matter 
power spectrum is measured directly from the cosmic shear 
data, upto a maximum k = SO/iMpc"^ . We find that relative 

to the non-marginalization case (where the power is assumed 
to be known exactly) there is a reduction in the ability of 
the cosmic shear survey to constrain dark energy parameters 
by 10%. This is equivalent to the maximum fc-range being 
reduced from k = 50/jMpc~^ to k ^ lOftMpc^^. 

We note that the information in the self-calibration 
regime comes from the geometric constraints on cosmolog- 
ical parameters that come from the cross-component to- 
mographic bin combinations (for a single bin, with auto- 
correlation information only the marginalized covariance, 
equation 11, —5- oo). Because the self-calibration re- 
covers the FoM to within 10% of the unmarginalized case, 
this suggests that the majority of the dark energy informa- 
tion from cosmic shear tomography, in the non-linear regime, 
comes from the geometric part of the lensing kernel. 



3.2 External Prior 

In this section we will place requirements on the functional 
variance from an external prior on the non-linear power spec- 
trum so that the expected dark energy cosmological con- 
straints remain unaffected. We parameterise the functional 
variance using 



al{k, 



P'^{k,z) 



■ Oo + aik y k> 0.5/iMpc 



(15) 



where we will investigate a constant uncertainty and one 
that scales linearly with the fc-mode, here k = £/r in equa- 
tion (14). This is an extension of the type of path integral 
marginalization used in Kitching & Taylor (2010) where only 
constant functional variances were considered, here we con- 
sider varying functional variance. Note that we assume that 
the variance is not a function of redhift, and that ai has 
units of /i~^Mpc. Figure 3 shows an example of the type of 
functional bounds that we will consider (we show the dimen- 
sionless power A^{k) = P{k)k^/2-K^ so that the non-linear 
regime is marked by A^(A;) > 1), throughout we assume that 
below k — 0.5/iMpc~^ the power spectrum is known exactly 
(zero functional variance) and we take an upper limit in 
wavenumber of A; < fcmax = 50/iMpc~^. We use the Eisen- 
stein & Hu (1999) linear power spectrum and the Smith et 
al. (2003) non-linear correction as the fiducial function. 

In Figure 4 we show how the FoM for a tomographic 
cosmic shear, for the survey outlined in Section 2, changes 
as a function of the functional variance on the non-linear 
power spectrum from an external prior. We find that if the 
functional variance can be constrained to < 1% then the 
FoM degradation is at most 1% relative to the case where 
no marginalization is performed. We provide a simple fit- 
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» 1— Sigma functional bound, aQ=0.5, aj = 0.0_ 

^ 1 0^ - 1— sigma functional bound, aQ = 0,5, a^ = G.O'o 
^ 1— sigma functional bound, ao=0.1, aj = 0.1 




0,001 0,010 0.100 1,000 10,000 100. OOC 
k [hMpcr'] 

Figure 3. The dark matter power spectrum at a redshift of zero, 
as a function of scale. We sliow 3 examples of the type of 1-sigma 
functional variances that can be considered using equation 15. 
Below k = 0.5/iMpc~^ we assume that the power spectrum is 
known exactly, which extends slightly into the non-linear regime. 
The horizontal dashed line marks A^(fc) = 1, which is where there 
is approximately a 1% deviation from linearity. For the calcua- 
tions in Section 3 that actual values of ao and ai are an order of 
magnitude smaller than those shown here, which are are illustra- 
tive purposes. 

ting formula for the dependency of the FoM on the constant 
contribution to the functional variance 

FoM = 224[1 - {Saof], (16) 
which is valid for ao <0.1. 



4 CONCLUSION 

To conclude we find that the non-linear and baryon- 
dominated part of the matter power spectrum contains, 
above wavenumbers of fc > 0.5/iMpc~^, half of the infor- 
mation content on dark energy parameters, parameterised 
through the Figure of Merit. However the lack of knowl- 
edge about this regime, and the complex simulated mod- 
elling needed to correctly constrain its behaviour as a func- 
tion of environment, scale and cosmology means that the 
uncertainty on the non-linear power spectrum must be cor- 
rectly accounted for in cosmic shear surveys. This article has 
some resonance with previous work, that parameterise the 
uncertainty in the non-linear power spectrum and marginal- 
ize over those parameters (Zhang et al., 2009; Rudd et al., 
2008; Zenter et al., 2008; Huterer et al., 2006; Jing et al., 
2000) or attempt to modify the data to minimise the effect 
(Huterer & White, 2005), however all these assume parame- 
terized models, that may not be able to reflect the real effect 
of baryons. 

We have derived likelihood expressions for the tomo- 
graphic cosmic shear power spectrum in the cases that the 
functional matter power spectrum is self-calibrated from the 
data itself, and in the case that an external prior on the func- 
tional variation of the matter power spectrum is available. 

We summarise our results in Figure 5. 

• With no external priors, a Euclid-like cosmic shear sur- 
vey, with a fcmax = 50/iMpc~^ could achieve a FoM ~ 220 



150 167 183 200 217 233 250 
FoM 

0.04 i^l^^l^l^B^I^S^^I 
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180 P I I I I ^ 

0.00 0.01 0.02 0.03 0.04 0.05 

Figure 4. The FoM, for a Euclid-like tomographic cosmic shear 
survey alone, as the functional variance of the external prior is 
varied, as parameterised by equation (15). The FoM degrades as 
the functional variance is increased. The upper panel shows how 
the FoM degrades as a functional of ao and ai, the lower panel 
show the dependency on ao keeping ai = 0. We provide a sim- 
ple fitting formula (blue, light gray, line) that relates a constant 
functional variance ao to the FoM. If the functional variance is 
< 0.01 (vertical dashed line) the FoM is degraded by less than 1%. 
The horizontal dashed line shows the FoM in the self calibration 
regime. 

from cosmic shear tomography alone. 

• If the non-linear matterpower spectrum is completely 
removed using a hard cut in fc-modes of fcmax ~ 0.5/iMpc~^ 
then the FoM is reduced by a factor of 50%. 

• In the functional self-calibration regime the cosmic 
shear survey can recover the FoM, with only a 10% reduc- 
tion in the FoM. 

• By including an informative prior, from simulations for 
example, the orginal FoM can be recovered if the func- 
tional variation of the non-linear matter power spectrum 
is known to ~ 1% to fc = 50/iMpc~^, or a physical scale of 
~ 120Kpc//i. 

Finally we note that in the self-calibration regime, the in- 
formation used to constrain the cosmological information, 
through the lensing kernel, has some similarities with the 
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Figure 5. Summary of the main results. If the non-Unear scales 
are removed using a hard cut in fc-modcs then the FoM can be 
reduced by 50%. In the functional self-calibration regime the re- 
duction is less severe with a relative reduction of 10%, finally if 
a 1% external prior on the functional variation of the non-linear 
power can be applied then the FoM is recovered. All fc-modes are 
in units of [ft.Mpc~^]. The parameters ao and ai refer the shape 
of the functional variance as a function of scale, given by equation 
(15). 

shear-ratio method, where cluster scale weak lensing is iso- 
lated. 

Constraining the non-hnear power spectrum to 1% func- 
tional accuracy down to 120kpc//i as a function of scale, 
redshift and cosmology is a significant theoretical and ob- 
servational challenge. On the theoretical side modeling the 
baryons on the scale of galaxy clusters is already a chal- 
lenge. Extending this to group and individual galajcy haloes 
will require a much deeper understanding of the baryonic 
processes on these scales. On the observational side, weak 
lensing itself, and galaxy-galaxy lensing, can provide much 
empirical information about the mass distribution which can 
be compared with stellar and gaseous components. These are 
challenges that must be realised if we are to fully exploit the 
potential of tomographic cosmic shear experiments. 
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